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Abstract 

We consider pairing in a three-component gas of degenerate fermions. In particular, we solve 
the finite temperature mean-field theory of an interacting gas for a system where both interac- 
tion strengths and fermion masses can be unequal. At zero temperature we find a a possibility 
of a quantum phase transition between states associated with pairing between different pairs of 
fermions. On the other hand, finite temperature behavior of the three-component system reveals 
some qualitative differences from the two-component gas: for a range of parameters it is possible 
to have two different critical temperatures. The lower one corresponds to a transition between dif- 
ferent pairing channels, while the higher one corresponds to the usual superfluid-normal transition. 
We discuss how these phase transitions could be observed in ultracold gases of fermionic atoms. 
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I. INTRODUCTION 



Recently, the possibility of fermion pairing in a two-component system with unequal 
Fermi surfaces has attracted a considerable amount of attention ^ . In an electronic system 
non-matched Fermi surfaces could be due to the magnetic field interacting with electron spins 
and in color superconductivity, due to unequal quark masses. The newly reahzed strongly 
interacting superfluid Fermi gases [1, 0, B S & H S offer a promising playground for 
the study of pairing and superfluidity, also with imbalanced Fermi energies. 

In ultra cold degenerate gases, different components are typically atoms in different in- 
ternal states of an atom, but different isotopes can also be considered. In ultra cold gases 
the mismatch of the Fermi surfaces can be due to having an unequal number of atoms in 
different states. After the first experimental studies using a mixture of two internal states 
of ^Li atoms [H, the properties of such systems, in a harmonic trapping potential, have 
been extensively discussed in the recent literature [3, 3 3, 3, 12) 3 13 2^ 21. 22. 23l|. 

In principle, many different types of fermions can be trapped in a same trap. In order 
to reach long lifetimes, the atoms should have favorable collisional properties, but current 
knowledge of these collisional properties is limited. However, mixtures of either Bose con- 
densed or fermionic atoms in a variety of different internal states have been experimentally 
demonstrated 2^, 2^ 2^ 27, 2^, and experimental effort towards realizing three-component 
Fermi gases has been initiated. For this reason we explore the novel possibilities in a three- 
component Fermi gas by generalizing the BCS (Bardeen-Cooper-Schrieffer) theory into three 
components. 

Related three-component Fermi system has been previously studied by Honerkamp and 
Hofstetter H ^ and more recently by Bedaque and D'Incao jsj]. The first reference 



employs a Hubbard lattice Hamiltonian and is focused on the rather special case when all 
fermions have the same mass and interact with different fermions (as well as with fermions 
of the same type) with a single interaction strength. The setting used in the latter reference 
includes the possibility of unequal interaction strengths between various components. Be- 
daque and D'Incao draw general qualitative conclusions on the symmetries of the possible 
zero temperature phases, but all fermion masses are assumed equal. 

In this paper we discuss pairing in a three-component system both at zero as well as 
non-zero temperatures. Furthermore, we allow for the possibility that the third fermion 
type has a different mass from the other two and since a priori there is little reason to 
expect identical scattering properties between different fermion components, we formulate 
our theory assuming different interaction strengths between the components. In ultracold 
gases, such a system could be realized by trapping, for instance, ^Li and ^^K using atoms in 
two internal states for one of theses and one internal state for the other In order to restrict 
the parameter space somewhat, we make one simplifying assumption in our Hamiltonian. 
Generally one expects three different interactions in a three-component system. One between 
the first and the second component, one between the second and the third component, and 
one between the third and the first component. In this paper we assume that the interaction 
between the third and the first component is weaker than the other two interactions and 
can be ignored next to the dominant contributions. This restriction is simply a matter of 
convenience as it restricts parameter space to manageable proportions, but does not affect 
the qualitative picture we find. 

We find that, within the BCS theory, the possible pairing always gives rise to just one 
order parameter of broken U{1) symmetry. Whether the pairing takes place between the first 
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and the second component (1 — 2 channel) or between the second and the third component 
(2 — 3 channel) depends on the strengths of the interactions, atomic masses, as well as 
on the differences between the chemical potentials of different components. The unpaired 
component constitutes a normal Fermi sea. At zero temperature there is a possibility of a 
quantum phase-transition from a phase with pairing in the 1 — 2 channel into a phase with 
pairing in the 2 — 3 channel as chemical potential differences are varied appropriately. This 
transition is of the first order. Furthermore, as the temperature increases it is possible to 
have a second order transition from a paired phase in the 2 — 3 channel into a paired phase 
in the 1 — 2 channel. At even higher temperature, there is another transition from the paired 
state into the normal state. 

This paper is organized as follows. In Sec. |n] we discuss the normal three-component 
Fermi gas and the possibility of phase-separation in such a system. In Sec. IIIII we present 
the BCS-style mean-field theory and discuss the quahtative behavior of the associated grand 
potential at zero temperature. In this section we present the energy landscape of the grand 
potential for several different scenarios. In Sec. lIVI we find the solutions to the gap-equations 
that correspond to the global minimum of the grand potential for fermions where the third 
component has a different mass from the other two components. We end with some con- 
cluding remarks in Sec. 



II. NORMAL THREE COMPONENT FERMI GAS 

We study a homogeneous three component Fermi gas whose components are either atoms 
in different internal states or atoms of different isotopes. For now we assume that an in- 
teraction between the first and the the third component is sufficiently small, so that it is 
justified to focus on only the interactions between the first and the second component and 
the second and the third component. The second quantized Hamiltonian for this system is 
therefore 

\cr=l,2,3 \ o- / 

+ 912 J rfrV^|(r)V^|(r)V^2(r)V^i(r) +^23 J c^r V^3(i")^2(r)^2(r)^3(r), 

where ipai'^) and '^/'^(r) are the usual field-operators which annihilate and create particles 
in state a. In addition, is the atomic mass, gi2 is the interaction strength between the 
first and the second component, and ^23 is the interaction strength between the second and 
the third component. In terms of the scattering lengths Qij these interaction strengths are 
given by gij = 27c?i^aij/ fj., where fi = mimj/{mi + mj) is the reduced mass of the scattering 
atoms. 

Let us first discuss a mixture of three normal fermion components. If the mixture of 
normal three-component gas phase-separates spatially, it is clear that there is no possibility of 
a superfluid where all three components are mixed together and the problem is that of a two- 
component Fermi gas surrounded by the normal third component. Since the temperatures 
of interests are very low compared to the Fermi energies of various components we can 
consider the zero-temperature limit, in which the free-energy density / of the mixture of 
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normal Fermi gases is given by 



2^2 ^5/3 ^5/3 ^5/3 

/ = ^(6vr')'/'(^ + ^ + ^) + ~gi2nin2 + ^23n2n3, (2) 
10 nil ^2 "^3 

where are the component densities. 

If both interaction strengths ^12, ^23 are positive then all three components separate. The 
reason for this is that then the positive definite free-energy Eq. (0) is minimized when the 
interaction terms give vanishing contributions. If one of the interactions is positive and the 
other is negative, the components with attractive interaction mix and the third component 
separates. Therefore the case when all three components can coexist spatially occurs only 
when both interactions are attractive. If one of the interactions is attractive then formally 
Eq. (j2)) has a global minimum at infinite density of the interacting components. This global 
minimum is not physically relevant since there is an energetic barrier separating this solution 
from the physically relevant regime of densities. In fact the height of this barrier increases as 
the interactions become weaker. From now on we will focus on the most interesting region 
where both interactions are fairly large and negative. 



III. GRAND POTENTIAL OF THE THREE COMPONENT FERMI GAS 

The BCS-theory involves approximating the Hamiltonian by the mean-field Hamiltonian 
density (in the k-space) 

HbCS = T7 E E ( 2 I ^'^kV^'^k + Ai2^i\kV^2,-k + ^l2V^2,-kV^l,k + 

k <x=l,2,3 V / 



lA |2 lA |2 

',t J,t I A* J,_ . . 1 ^121 1^231 

9l2 932 



+ ^23^3,k^2,-k + A23V'2,-kV'3,k 



where 



are the order parameters to be determined self-consistently. As a short hand notation we 
define the vector A = (A12, A23) characterizing the state of the system. 

Since the mean-field Hamiltonian is of second order in the operators it can be easily diag- 
onalized with a canonical transformation. In this way we can calculate the grand potential 
(or free energy) 

^](Al2, A23) = -kBT\og[Tr{exp{-(3HBcsm (4) 

where ks is the Boltzmann constant and f3 = l/fc^T. In the continuum limit the grand 
potential is ultraviolet divergent. This divergence is caused by the unphysical short distance 
behavior of the contact interaction and is removed, in the usual way, by subtracting the 
divergent contribution from the grand potential. We have done the diagonalization of the 
mean field Hamiltonian which can be done analytically, but in the general case formulas are 
inconveniently long and not very informative. In order to get an overview of the expected 
system behavior it is rather more instructive to focus on the behavior of the free-energy. 

We choose our units by assuming that the first and the second components are ^Li atoms 
and by using the ideal gas of the second component as a benchmark. The unit of energy 
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is then ep = h'^kp2/'^''^2, where kp2 = (671^^2)^/^ is the Fermi wave-vector. Furthermore, 
we use the unit of length l/kp^2- Using these units we define the dimensionless couphng 
strengths as gi2 = 2kF,2ai2/''T and g23 = (m2 + ■m3)kF,2a23/'m'2'iT. 

A. Grand potential when all fermions have the same mass 

The Fig. Q demonstrates the behavior of the free energy landscape as a function of the 
gaps A12 and A23 when all fermion components have the same mass. It is clear that when 
all chemical potentials are equal the free energy landscape is only a function of Af2 + ^^23 
and at zero temperature displays a clear minimum corresponding to the solution of the 
B CS gap equ ation. With the parameters used in the figure, the minimum is located at 
a/ + ^23 ~ 0.135, see the quarter sphere in Fig. [U^a). 

Any difference in the chemical potentials (or coupling strengths) breaks the above sym- 
metry and makes it energetically favorable to have another one of the phases A = (A12, 0) 
or A = (0,A23). When interactions are the same, which of these alternatives is chosen 
depends on the average chemical potential between the paired components. The channel 
corresponding to the higher average chemical potential has a lower free energy and is there- 
fore physically realized. At zero temperature, the densities of the paired components are 
equal, as is to be expected from the BCS-state. At higher temperatures the normal state 
A = (0, 0) eventually becomes the free energy minimum. The Figure ^ also demonstrates 
that generally the grand potential has several possible stationary points. 

B. Grand potential when the fermion masses differ 

The different components could in principle have also different masses. For concreteness 
we assume that the first and the second component are different internal states of ®Li atoms 
while the third component is a ^''K atom. This therefore corresponds to the mass ratio 
rrir = rrii/m^ = 0.15. In the Fig. |21we show some typical free energy landscapes in this case, 
when the chemical potential of the third component is varied. It is clear that the grand 
potential again has many different stationary points. 

For small values of the pairing is only possible in the 1 — 2 channel. However, the 
Fermi surfaces are matched when /is = and this is reflected as a possibility of pairing in 
the 2 — 3 channel around this value. The reason why pairing in this channel is favored for 
matched Fermi surfaces even when coupling strengths are equal, is due to the higher density 
of states for atoms of higher mass. This higher density of states translates into reduction in 
energy. Naturally if (723 is reduced sufficiently we enter a parameter region where pairing in 
the 2 — 3 channel will never take place. For the parameters used in the figure this happens 
when |5f23| < 0.33. 

Again we find that at zero temperature the densities of the paired components are always 
equal. It is also interesting to observe that by increasing the temperature, one can induce 
a transition from the A = (0, A23) phase into the A = (Ai2,0) phase before entering a 
normal state. 

In summary, two transitions are visible: a zero temperature quantum phase transition 
between the two pairing channels (transition from Fig. |2fa) to Fig. |2fb), and from Fig. |2fb) 
to Fig. |2fc)), and a finite temperature second order transition between the pairing channels 
(transition from Fig. IH^b) to Fig. I21d)). 
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FIG. 1: (Color online) The free energy landscape as a function of A12 and A23 when all three 
components have the same mass. The dash-dotted line in (a) and the filled circles in (b)-(d) 
show the location of the global minimum. We used equal coupling strengths gi2 = (723 = —0.50 
and the figures (a)-(c) were calculated at zero temperature while the figure (d) was calculated at 
ksT/ep = 0.08 which is above the critical temperature. The chemical potentials we such that in 
the figure (a) /Ui = /U2 = /X3 = 1, in (b) ^ui = 1.01, fi2 = 1, /^3 = 0.99, in (c) /ii = 0.99, ;U2 = 1, 
/i3 = 1.01, and finally in (d) jJ-i = H2 = IJ-'S = 1- 



IV. SOLUTIONS OF THE GAP-EQUATIONS OF THE THREE COMPONENT 
FERMI GAS 

For a given set of chemical potentials the self-consistent BCS-solution is found at the 
minimum of the grand potential. At the extremum value of the grand potential, the gap- 
equations 

gr](Al2,A23) 

and 

gr](Al2,A23) 
9A23 



(5) 
(6) 
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FIG. 2: (Color online) The free energy landscape as a function of A12 and A23 with the mass 
ratio mi/m3 = 0.15 and = ^2 = l-O. The filled circle shows a location of the global minimum. 
The coupling strengths were gi2 = 923 = —0.5 and the figures (a)-(c) were calculated at zero 
temperature while the figure (d) was calculated at ksT/ep = 0.07. In figure (a) fi^ = 0.07, in 
(b) Us = 0.15, in (c) = 0.25, and in (d) //s = 0.15. (In the figure (d) we used a logarithmic 
scale to calculate the contour plot in order to enhance the relevant features.) The results show 
two types of transitions: a zero temperature quantum phase transition between the two pairing 
channels (from (a) to (b), and from (b) to (c)), and a finite temperature second order transition 
between the pairing channels (from (b) to (d)). 

are satisfied. These equations are satisfied at all extremal points of the grand potential, 
including local minima, local maxima, and saddle points. As we found out in the earlier 
section, the grand potential can have several extremal points so care must be taken in 
ensuring that numerics converges to the physically relevant global minimum of the grand 
potential rather than to a local stationary point. 

We have found the global minima of the free energy by solving the gap equations under 
many different circumstances. Fig. |21 shows typical results for the gaps as a function of 
temperature and fi^. At zero temperature there exists a first order quantum phase transition, 
as the chemical potential is varied, from the A = (A12, 0) phase into the A = (0, A23) phase. 
We also find that with increasing temperature there can be a second order transition from 



7 



(a) (b) (c) 




FIG. 3: The gap parameters as a function of temperature T and of the thh'd component chemical 
potential 113. We used the mass ratio rrii/m^ = 0.15, /Ui = /i2 = 1, and the coupling strengths 
were 512 = 523 = —0.5. The figure (c) shows A = Af 2 + and demonstrates clearly the sharp 
change in the order parameter at zero temperature (the quantum phase transition) as well as the 
smoother transition from the A23 superfluid to A12 superfluid at a finite temperature. 

the A = (0, A23) phase into the A = (A12, 0) phase, which then makes a transition into the 
normal state at a higher temperature. However, whether this sequence of transitions takes 
place or not, does depend on the strength of the interactions. For weaker interactions it 
is possible that the A = (0, A23) phase is not completely surrounded by the A = (Ai2,0) 
phase as it is in Fig. El 

This overall behaviour can be understood in the following way: The pairing between 
atoms of different masses happens around the point where their Fermi momenta are nearly 
equal (following from the condition 7x3/ rrii/m^). This pairing is favoured over the pairing 
with same mass components due to the higher density of states for the more massive compo- 
nent. At finite temperature, the Fermi surfaces of both pairing components will be smeared 
out due to temperature - however, for components with different masses, this smoothing 
of the Fermi edge will happen differently (the higher mass component Fermi edge spreads 
more). Therefore, with increasing temperature, it becomes increasingly difficult to find a 
pairing partner with the same magnitude (but opposite sign) momentum, and also the ad- 
vantage given by the higher density of states of the more massive component becomes less 
significant. This makes the pairing between the same mass components potentially more 
favourable at high temperatures. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have considered the mean-field theory of an interacting three-component 
Fermi gas. We found that in the highly symmetric case when all fermion masses as well as 
interactions strengths are equal the free energy is only a function of A^2 + A23- However, any 
deviation from the most symmetric situation, whether it is due to different atomic masses, 
different interaction strengths, or different chemical potentials, breaks the above symmetry 
of the free energy and makes one of the paired states (A12 7^ 0, A23 = 0), (A12 = 0, A23 7^ 0) 
or the normal-state (A12 = 0, A23 = 0) energetically favored. Magnitudes of non-zero gaps 
as well as the point where the transition to the normal-state occurs depend on parameters, 
but we have solved the global minima of the free energy for fairly representative range of 
parameters. 

When using a grand canonical ensemble the average densities are quantities derived from 
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the grand potential through = —dfl/dfii. In this work we have explored the system 
behavior as a function of the chemical potentials by minimizing the grand potential. It is 
therefore useful to elaborate what kind of density ratios different sets of chemical potentials 
actually correspond to. For the case investigated in the earlier section with a mass ratio 
mi/m^ = 0.15 at zero temperature, we find at /is = 0.05 pairing in the 1 — 2 channel and 
density ratios of n^/rii = n^/n2 ~ 0.26. As the third component chemical potential is 
increased to /xa = 0.15, we find pairing in the 2 — 3 channel with density ratios n^jnx = 
n^/ni ~ 0.72. Finally, when /i3 = 0.25, the pairing is reverted back to the 1 — 2 channel 
and the mass ratios are n^/ni = n^/n2 ~ 2.9. Therefore, this transition could be realized 
simply by varying the density of one component and keeping the other two fixed. Note that 
transitions could also be induced by changing the interaction strengths, i.e. the scattering 
lenghts, between the three components as discussed in [31||. However, at the present the 
method used for tuning the scattering lengths, i.e. the use of Feshbach resonances, does 
not allow to tune the different scattering lengths independently whereas the densities can 
be varied completely independent of each other. The change in the order ga p re lated to 
the transitions could be observed using RF-spectroscopy of the pairing gap jsl. l32l 33l | and 
corresponding condensate fraction by the pair projection method j^b iSj, and superfluidity 
could be directly demonstrated by the creation of vortices |l3| - all these can be studied 
selectively for any choice for the pairing between the three components. 

In this paper we have for convenience limited ourselves to the special case where inter- 
action between the first and the third component is small enough to be neglected. This 
requirement can be relaxed easily and one could solve the most general case with unequal 
masses as well as with unequal interactions. We have checked that the presence of the omit- 
ted interaction between the first and the third component does not change the qualitative 
picture presented in this paper. For some parameter values the new gap-function A13 can 
become non-zero, but otherwise the qualitative physics remains similar to the cases discussed 
in this paper. 



Three component Fermi system has been studied theoretical |29l l30l l31l ? ]. The first 
reference employs a Hubbard lattice Hamiltonian and is focused on the rather special case 
when all fermions have the same mass and interact with different fermions (as well as with 
fermions of the same type) with a single interaction strength. The setting used in the 
latter reference includes the possibility of unequal interaction strengths between various 
components. Bedaque and D'Incao draw general qualitative conclusions on the symmetries 
of the possible zero temperature phases, but all fermion masses are assumed equal. 
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